Gluon propagators and center vortices in gluon plasma 
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We study electric and magnetic components of the gluon propagators in quark-gluon plasma in 

m ; 

terms of center vortices by using a quenched simulation of SU (2) lattice theory. In the Landau 

■ 

\ gauge, the magnetic components of the propagators are strongly affected in the infrared region by 



o 



removal of the center vortices, while the electric components are almost unchanged by this pro- 
cedure. In the Coulomb gauge, the time-time correlators, including an instantaneous interaction, 



also have an essential contribution from the center vortices. As a result, one finds that magnetic 



degrees of freedom in the infrared region couple strongly to the center vortices in the deconfinement 
phase. 
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I. INTRODUCTION 



The Relativistic Heavy- Ion Collider (RHIC) at Brookhaven National Laboratory produces 
a new state of matter which may exceed the critical temperature T c of the phase transition 
from the hadron phase to the quark-gluon plasma (QGP) phase. Many phenomenological 
studies and lattice computations suggest that the QGP is a strongly interacting plasma jlj], 
for which we cannot apply the early arguments based on the perturbative approach with a 
small coupling constant. Furthermore, the recent Pb-Pb heavy- ion collision experiment at 
LHC has created QGP matter at even higher temperatures: This shows us an obvious jet- 
quenching event 2j and a larger elliptic flow [3| compared to the RHIC's Au-Au collision. 
Therefore, it is indispensable for us to explore the mechanism which drives such strong 
interactions using a nonperturbative first-principle approach in lattice simulations. 

One of the most important ideas to describe a strongly interacting QGP (sQGP) is to fo- 
cus on an infrared singularity arising from magnetic degrees of freedom {4, sj]. The magnetic 
component of the gluon propagator is fully inaccessible by the perturbative calculation, but 
its infrared divergence may cause an emergence of a nonperturbative magnetic mass that 
olays a cutoff role and can cure thermal QCD in the infrared region. The lattice simulations 
^-8] prove that the magnetic gluons have a nonvanishing mass at finite temperature. Fur- 
thermore, it is well known that a spatial Euclidean Wilson loop (which is not extended to 
the temporal dimension) bears a confining potential above T c 9l4l7j]. while the correlators of 
a Polyakov line - wrapped in the temporal direction - give a nonconfining screened potential 
of the Debye type, with a finite electric mass ~ g{T)T 



17 



20|. 
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In addition, from the viewpoint of the Gribov-Zwanziger (GZ) confinement scenario 
23|, a color-Coulomb instantaneous interaction between a quark and an antiquark provides 
- even in the nonconfined QGP phase - a confining potential which rises linearly as the 



function of the quark-antiquark separation 2J- 



271 ] . As a consequence, the thermal string 



tensions obtained from the spatial- Wilson and the color-Coulomb potentials are nonzero. 
They depend on the temperature and obey a magnetic scaling law [~ g 2 (T)T]. Extending 
this line of considerations, Zwanziger has approximately reconstructed the equation of state 



32|. 



of QGP using the Gribov-type dispersion relation for the massive gluons 

These interesting aspects of the non-Abelian gauge theory may be related to center (mag- 
netic) vortices - i.e., to the topological defects associated with the nontrivial homotopy group 



7Ti[SU(N) /Z(iV)] ~ Z(JV) - which are responsible for certain nonperturbative phenomena 
of QCD. One can identify the center vortices on the lattice using a numerical technique 33 1 
and also remove these vortices from the original gauge fields 35|. It turns out that the re- 
moval of the center vortices destroys the color confinement property and restores the chiral 
symmetry. Moreover j_the lattice center- vortex density exhibits a scaling consistent with the 



asymptotic freedom 



In terms of the vortex degrees of freedom, the QCD deconfinement phase transition can be 
considered as a depercolation transition of the vortex lines in the direction of the Euclidean 

n 

time [37|. As a result, we can naturally understand the survival of the spatial confinement 
above T c because the center vortices remain intact in the spatial space. Moreover, a typical 
center-vortex configuration is located at the Gribov horizon in the gauge space. Thus, 
the removal of the center vortices results in the dilution of the lowest eigenvalues of the 
Faddeev-Popov operator. These eigenvalues - according to the GZ confinement scenario - 
cause confinement of color |25, 



34]. 



Recently, three of us have argued that the center-vortex mechanism is also important in 



theory because the center vortices carry information about 



the hot phase of the Yang-Mills 
the magnetic degrees of freedom [38|, |39j. The center vortices are related to Abelian magnetic 
monopoles, and the later are expected to explain some of the interesting properties of the 
quark- gluon plasma as well [40 ]. 



In this paper we study a connection between the center vortices and the infrared properties 
of the gluon propagators at finite temperature. To this end we study the behavior of the 
electric and magnetic components of the gluon propagators by removing the vortices from 
the original gauge configurations and comparing the result with the original one. We use 
the quenched SU (2) lattice simulations in the Landau and Coulomb gauges. In Sec. II we 
define gluon propagators on the lattice. In Sec. Ill a numerical technique used to make a 
center projection is summarized. Our numerical results are presented in the Sec. IV, while 
the last section is devoted to the summary of this work. 
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II. GLUON PROPAGATORS 



In this study, we work in the 377(2) lattice gauge theory. The gauge potential A is 
expressed via the SU(2) matrix link variable C/ /J ( follows: 

A li {x,t) = ^^2 r Tra a U^x,t), W 

a 

where a a are the Pauli matrices. The correlation functions of the gauge fields (pQ) in mo- 
mentum space are 

D»M,t) = ^J2^^)My,t"))e iq{x - y \ (2) 

x,y 

where V(= N x N y N z ) is the three-dimensional volume and t — t' — t" is the Euclidean time 
difference. 

In Landau-gauge fixing we study the static correlators of gluon fields with g = 0: 

D^q, Qo = 0) = ^Yl D /-(<7, *) , (3) 

* t 

where N t is the lattice size in the Euclidean temporal direction. In the Coulomb gauge, it 
is more appropriate to investigate an equal-time gluon propagator in the following form: 

D%(q) = ^J2(M^t)My,t))e^\ (4) 

This propagator corresponds to the one in Eq. (J2]) at t' = t" . Note that there is no go 
dependence in Eq. ([2]) and that the go = term is removed from the sum. An equal-time 
propagator reads D(q) = l/(2u(q)) where uj = \J q 2 + m 2 is the dispersion relation. 

In the finite-temperature system, the electric and magnetic gluons have different effects 
due to breaking of the Euclidean Lorentz invariance. One can define the spatially transverse 
(Pt) and spatially longitudinal (Pl) projection operators as follows: 

pOO _ pOi _ piO pij _ rij 1 ^ (r\ 
i ^ — T — J T ' T — 2 ' \ / 

p»u = 5 ^ _<£?__ pjw ^ (g) 

with the properties 

(P T ) 2 = P T , (P L ) 2 = P L , P t Pl = 0. (7) 
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Both spatially transverse and spatially longitudinal projectors correspond to the transverse 
states in momentum space: 

%k v = %n v = o . (8) 

Using these relations, the gluon propagators at finite temperature in a Landau-type gauge 
can be separated into two independent parts: 

The electric component of the gluon propagator is given by the spatially longitudinal pro- 
jection, D e = D o and the electric mass is given by F((q,q ) = 0) = m 2 ~ (g(T)T) 2 . The 
spatially transverse projection gives us the magnetic propagator D m = Da. The magnetic 
mass is expected to be G((q,qo) = 0) = rr? m ~ (g 2 (T)T) 2 , where g(T) is a running QCD 
coupling defined at the scale of temperature T. 

III. MAXIMAL CENTER PROJECTION 

We employ a direct maximal center projection (MCP) [3^] in order to identify the center 
vortices on the lattice. The corresponding gauge is defined by the condition 

Maximize R = ^ Tr [U^(x)f . (10) 

The center gauge field, 

Z M (x) = sgnTr[^(x)]GZ 2 , (11) 

allows us to identify the center vortices. If the center plaquette is not equal to a trivial 
element (unity) then a center vortex goes through this plaquette. 



In order to remove the center vortices from the gauge-field ensemble, we follow Ref. 35 1 
by multiplying the original field by the center-projected field Z^. 

[/» = ZMx) , (12) 

so that the new links V ' correspond to vortex-free ensembles. 

It is confirmed by lattice simulations that the confinement and chiral symmetry breaking 
are both lost after the removal of the center vortices 33|, |35j. We would also like to note that 



the effect of the vortex removal on chiral symmetry breaking (%SB) depends on the choice of 
the lattice quark action 4lM43|. and thus chiral symmetry breaking should be treated with 
care. In our paper we address the problem of color confinement. 



IV. NUMERICAL RESULTS 



A. Lattice setup 



We carry out quenched SU{2) lattice simulations by generating gauge configurations. The 
convergence criterion of the MCP technique is set as 10 -16 , and the precision of an iterative 
gauge fixing algorithm |44j is 10 -8 . We use the single geometry of the lattice, 24 3 x 4, and 
various temperatures T/T c ~ 1.4, 3.0, and 6.0 [T c ~ 305 MeV for N t = 4 for the SU(2) 



gauge group] for (3 = 2.40, 2.64, and 2.88, respectively [39 



451 ] . We used approximately 20 



to 30 lattice configurations collected every 100 sweep steps. 



B. Thermal gluon propagators 

In the left plot of Fig. [1] we show the gluon propagators calculated in the Landau gauge 
at T/T c = 1.40. The removal of center vortices visibly affects the infrared behavior of the 
magnetic and electric gluons in the infrared region. However, the effect of the vortex removal 
is much more pronounced for the magnetic degrees of freedom compared to the effect on the 
electric correlators. 

The effect of the vortex removal on the electric component of the gluon propagator dimin- 
ishes with an increase of temperature, as one can see from the plots of Fig. (TJ corresponding 
to the higher temperatures T/T c = 3.0 and 6.0. However, the magnetic propagators are 
affected drastically by the center vortices in the infrared region for all studied temperatures. 

The Coulomb gauge gluon propagators are plotted in Fig. |2j The magnetic propagator is 
affected by the removal of the center vortices in the deconfinement regions, being consistent 
with that of the Landau-gauge case as we discussed in the previous paragraph. Contrary to 
this fact, the electric parts in the Coulomb gauge are influenced by the magnetic vortices. 
This tendency remains for T/T c = 1 ~ 6; thus, it seems that there is an inconsistency 
between the two gauges. 

However, we have to mention the fact that in the Coulomb gauge the temporal-gauge 
correlator is dominated by an infrared singularity arising from spatial (magnetic) compo- 



nents 
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2l-l30l] . The time-time correlator with Coulomb-gauge fixing can be decomposed into 
two parts: 

D 00 (x,t) = V c (x)5(t) + P(x,t), (13) 
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FIG. 1. Electric and magnetic gluon propagators in the Landau gauge for T/T c = 1.4, 3.0, and 6.0 
as a function of the spatial momentum. The circle symbols represent numerical results obtained 
with the original lattice configurations, while the square symbols correspond to the vortex-removed 
configurations. 



where V c is an instantaneous potential, which is responsible for the color confinement and 
P corresponds to the vacuum (retarded) polarization term. In this theory, V c is related to 
the Green's function M~ l of the Faddeev-Popov ghost, 

V c (x - y)5 ab = ((M-\-d*)M-Xy) , (14) 

which does not explicitly depend on the Euclidean time (temperature), and thus, this quan- 
tity has no effect on screening. In contrast, P is a function of time and it may contribute 
to the screening. Indeed, in the deconfinement phase, this term provides a screened quark 
potential with finite electric mass. The screening can be observed by the investigation of 



a Polyakov line correlator |8|. The Polyakov line correlator with Landau-gauge fixing gives 
the color-screened potential as well (3, Q . 

In the confinement region, V c is a linearly rising potential. Moreover, even above the 



critical temperature T c the potential V c is a confining potential 



251 ] . Its thermal color- 



Coulomb string tension depends on temperature. The temperature dependence is consistent 



with the magnetic scaling g 2 (T)T 27|. The remnant confinement property corresponds 
to the nonvanishing spatial string tension. Consequently, it is now obvious that the time 
correlator in the Coulomb gauge is also a magneticlike quantity. 

Conversely, the covariant-type Landau gauge may not plainly separate the longitudinal 
and transverse modes. Actually, it is more difficult to observe a confining property of gluons 



even in the confinement region, compared to the case of the Coulomb gauge 



311 ] . However, as 



seen in Fig. [TJ the physical magnetic gluon is definitely affected by the vortex removal, while 
the corresponding variation of the electric gluon (excluding the instantaneous interaction 
which is singular in the infrared region) is very small. 

It is natural that the confining behavior of the thermal gluon propagators has differ- 
ent forms for different gauge fixings. Nevertheless, it is very important to stress that the 
calculations in both gauges give us the same conclusion that the relevant elements to the 
magnetic degrees of freedom are strongly coupled to center vortices after the deconfining 
phase transition. 

In addition, our result means that the Gribov-Zwanziger confinement scenario survives 
above T c . According to this theory, the spatial correlator experiences the suppression effects 
in the confinement phase, and the temporal correlator diverges in the infrared limit. We see 
that a similar behavior is seen in our numerical data at finite temperature. Furthermore, our 
observation can also be derived from the fact that in the QGP phase the vortex configurations 
belong to the Gribov horizon 25]. 
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FIG. 2. Electric (upper figures) and magnetic (lower figures) gluon propagators in the Coulomb 
gauge as a function of the spatial momentum at T/T c = 1.4, 3.0 and 6.0. The open symbols 
represent the numerical results obtained with the use of the original lattice configurations, while 
the filled symbols correspond to the vortex-removed configurations. 
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C. Check for numerical ambiguities 



It is well known that the MCP has numerical ambiguities (the Gribov copies). In order 
to check the stability of our results against this ambiguity, we used random-gauge trans- 
formations (RGT) applied to the Monte-Carlo updated gauge-field configurations before 
performing the MCP. Although the global maximum of the gauge fixing functional ( TlO|) 
cannot be determined with an ideal accuracy, neither the electric nor the magnetic gluon 
propagator significantly depends on this algorithm, as is shown in Fig. [3j 



Q 




FIG. 3. The dependence of the electric (left panel) and magnetic (right panel) gluon propagators 
on the Gribov copy ambiguities. The lattice size is 24 3 x 4 and T/T c = 1.4. The RGT with different 
random seeds is applied to each lattice configuration a varying number of times. 



V. SUMMARY 

We have studied the effects of the center (magnetic) vortices on electric and magnetic 
components of the gluon propagators in the QGP phase using SU(2) lattice simulations in 
the Landau and Coulomb gauges. We find that the gluon dynamics in the infrared region 
strongly couple to the magnetic vortices. Thus, in the deconfinement phase the magnetic 
vortex degrees of freedom should be treated nonperturbatively. 

At high temperatures the removal of the magnetic vortices reduces drastically the mag- 
netic gluon propagators in the infrared region in both gauges. The effect is similar to the 

9 



suppression of the infrared gluon propagators in the confinement phase, observed first in 



Ref. 



46j |. The electric propagators are almost unaffected by this procedure in the Landau 



gauge while in the Coulomb gauge - in agreement with the GZ mechanism - the electric 
gluon propagator is suppressed by the removal of the center vortices. 
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